Assessing the performance of Moghani crossbred lambs derived from different mating systems with Texel and Booroola sheep.

In our ongoing project, which focuses on the introgression of Booroola/FecB gene and the myostatin (MSTN) gene into purebred Moghani sheep, we assessed the performance of second-generation Moghani crossbreds such as second crossbreds (F2) and initial backcross generation (BC1). These crossbreds were generated through different mating systems, including in-breeding, outcrossing, first paternal backcrossing (PBC1), and first maternal backcrossing (MBC1). Notably, F2 strains exhibited lean tail, woolly fleece and a higher percentage of white coat color compared to BC1. The impact of mating systems and birth types on pre-weaning survival rates was found to be statistically significant (P < 0.0001), with singleton offspring resulting from paternal backcross showing a particularly substantial effect. The F2 crossbred lambs carrying the Booroola gene did not show a statistically significant difference in survivability compared to those carrying the MSTN gene, implying the Booroola prolificacy gene had no significant impact on survival outcomes. However, the occurrence of multiple births had a significant negative impact on lamb survival (P < 0.0001). The PBC1 sheep strains, specifically Texel Tamlet ram strains carrying the MSTN mutation, exhibited superior growth rates compared to others (P < 0.05). Interestingly, the MSTN mutation in the homozygous variant genotype significantly impacts growth rate before weaning compared to other genotypes and pure Moghani sheep (P < 0.05). In conclusion, this study objectively underscores the pivotal role of genetic factors, specifically through strategic mating systems like paternal backcrossing, in enhancing desired traits and growth rates in Moghani sheep, thereby contributing valuable insights to the field of sheep breeding programs.

Searching for homozygous haplotype deficiency in Manech Tête Rousse dairy sheep revealed a nonsense variant in the MMUT gene affecting newborn lamb viability.

BACKGROUND: Recessive deleterious variants are known to segregate in livestock populations, as in humans, and some may be lethal in the homozygous state. RESULTS: We used phased 50 k single nucleotide polymorphism (SNP) genotypes and pedigree data to scan the genome of 6845 Manech Tête Rousse dairy sheep to search for deficiency in homozygous haplotypes (DHH). Five Manech Tête Rousse deficient homozygous haplotypes (MTRDHH1 to 5) were identified, with a homozygous deficiency ranging from 84 to 100%. These haplotypes are located on Ovis aries chromosome (OAR)1 (MTRDHH2 and 3), OAR10 (MTRDHH4), OAR13 (MTRDHH5), and OAR20 (MTRDHH1), and have carrier frequencies ranging from 7.8 to 16.6%. When comparing at-risk matings between DHH carriers to safe matings between non-carriers, two DHH (MTRDHH1 and 2) were linked with decreased insemination success and/or increased stillbirth incidence. We investigated the MTRDHH1 haplotype, which substantially increased stillbirth rate, and identified a single nucleotide variant (SNV) inducing a premature stop codon (p.Gln409*) in the methylmalonyl-CoA mutase (MMUT) gene by using a whole-genome sequencing approach. We generated homozygous lambs for the MMUT mutation by at-risk mating between heterozygous carriers, and most of them died within the first 24 h after birth without any obvious clinical symptoms. Reverse transcriptase-qPCR and western blotting on post-mortem liver and kidney biological samples showed a decreased expression of MMUT mRNA in the liver and absence of a full-length MMUT protein in the mutant homozygous lambs. CONCLUSIONS: We identified five homozygous deficient haplotypes that are likely to harbor five independent deleterious recessive variants in sheep. One of these was detected in the MMUT gene, which is associated with lamb lethality in the homozygous state. A specific management of these haplotypes/variants in the MTR dairy sheep selection program would help enhance the overall fertility and lamb survival.

Comparison of the growth performance between pure Moghani sheep and crosses with Texel or Booroola sheep carrying major genes contributing to muscularity and prolificacy.

This study aimed to compare growth performance between Moghani sheep and crossbred lambs resulting from crossbreeding between Moghani pure breed ewes and the lines of rams e.g., Texel Tamlet, Texel Dalzell, Booroola Merino, and Booroola Romney. The first visible phenotypic characteristic was the presence of lean tail in all F1 crossbred lambs, whereas Moghani pure sheep is a well-known large fat-tailed breed. Moreover, the first generation of backcross (BC1) lambs from mating four types of F1 crossbred rams with Moghani pure ewes revealed lean-tailed to short fat-tailed. Comparative results showed that the F1 crossbred lambs had significantly (p < 0.0001) greater birth weight (BW) than the Moghani pure breed lambs. Despite no significant differences observed between Moghani pure breed sheep and its F1 crossbred lambs for body weight at pre-weaning, but F1 crossbred lambs achieved significantly (p < 0.0001) greater body weight after weaning compared to Moghani sheep. The growth performance of BC1 lambs was outperformed than F1 crossbred and Moghani sheep. These results encourage the continuation of the Moghani sheep crossbreeding programs to improve overall lamb growth, particularly post-weaning and to benefit from a better reproductive efficiency by elimination or reduction of the fat tail.

At the crossroads of fertility and metabolism: the importance of AMPK-dependent signaling in female infertility associated with hyperandrogenism.

STUDY QUESTION: What biological processes are linked to the signaling of the energy sensor 5'-AMP-activated protein kinase (AMPK) in mouse and human granulosa cells (GCs)? SUMMARY ANSWER: The lack of α1AMPK in GCs impacted cell cycle, adhesion, lipid metabolism and induced a hyperandrogenic response. WHAT IS KNOWN ALREADY: AMPK is expressed in the ovarian follicle, and its activation by pharmacological medications, such as metformin, inhibits the production of steroids. Polycystic ovary syndrome (PCOS) is responsible for infertility in approximately 5-20% of women of childbearing age and possible treatments include reducing body weight, improving lifestyle and the administration of a combination of drugs to improve insulin resistance, such as metformin. STUDY DESIGN, SIZE, DURATION: AMPK signaling was evaluated by analyzing differential gene expression in immortalized human granulosa cells (KGNs) with and without silencing α1AMPK using CRISPR/Cas9. In vivo studies included the use of a α1AMPK knock-out mouse model to evaluate the role of α1AMPK in folliculogenesis and fertility. Expression of α1AMPK was evaluated in primary human granulosa-luteal cells retrieved from women undergoing IVF with and without a lean PCOS phenotype (i.e. BMI: 18-25 kg/m2). PARTICIPANTS/MATERIALS, SETTING, METHODS: α1AMPK was disrupted in KGN cells and a transgenic mouse model. Cell viability, proliferation and metabolism were evaluated. Androgen production was evaluated by analyzing protein levels of relevant enzymes in the steroid pathway by western blots, and steroid levels obtained from in vitro and in vivo models by mass spectrometry. Differential gene expression in human GC was obtained by RNA sequencing. Analysis of in vivo murine folliculogenesis was performed by histology and immunochemistry, including evaluation of the anti-Müllerian hormone (AMH) marker. The α1AMPK gene expression was evaluated by quantitative RT-PCR in primary GCs obtained from women with the lean PCOS phenotype (n = 8) and without PCOS (n = 9). MAIN RESULTS AND THE ROLE OF CHANCE: Silencing of α1AMPK in KGN increased cell proliferation (P < 0.05 versus control, n = 4), promoted the use of fatty acids over glucose, and induced a hyperandrogenic response resulting from upregulation of two of the enzymes involved in steroid production, namely 3ß-hydroxysteroid dehydrogenase (3ßHSD) and P450 side-chain cleavage enzyme (P450scc) (P < 0.05, n = 3). Female mice deficient in α1AMPK had a 30% decrease in their ovulation rate (P < 0.05, n = 7) and litter size, a hyperandrogenic response (P < 0.05, n = 7) with higher levels of 3ßHSD and p450scc levels in the ovaries, and an increase in the population of antral follicles (P < 0.01, n = 10) compared to controls. Primary GCs from lean women with PCOS had lower α1AMPK mRNA expression levels than the control group (P < 0.05, n = 8-9). LARGE SCALE DATA: The FastQ files and metadata were submitted to the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB46048. LIMITATIONS, REASONS FOR CAUTION: The human KGN is a not fully differentiated, transformed cell line. As such, to confirm the role of AMPK in GC and the PCOS phenotype, this model was compared to two others: an α1AMPK transgenic mouse model and primary differentiated granulosa-lutein cells from non-obese women undergoing IVF (with and without PCOS). A clear limitation is the small number of patients with PCOS utilized in this study and that the collection of human GCs was performed after hormonal stimulation. WIDER IMPLICATIONS OF THE FINDINGS: Our results reveal that AMPK is directly involved in steroid production in human GCs. In addition, AMPK signaling was associated with other processes frequently reported as dysfunctional in PCOS models, such as cell adhesion, lipid metabolism and inflammation. Silencing of α1AMPK in KGN promoted folliculogenesis, with increases in AMH. Evaluating the expression of the α1AMPK subunit could be considered as a marker of interest in infertility cases related to hormonal imbalances and metabolic disorders, including PCOS. STUDY FUNDING/COMPETING INTEREST(S): This study was financially supported by the Institut National de la Recherche Agronomique (INRA) and the national programme « FERTiNERGY ¼ funded by the French National Research Agency (ANR). The authors report no intellectual or financial conflicts of interest related to this work. R.K. is identified as personnel of the International Agency for Research on Cancer/World Health Organization. R.K. alone is responsible for the views expressed in this article and she does not necessarily represent the decisions, policy or views of the International Agency for Research on Cancer/World Health Organization. TRIAL REGISTRATION NUMBER: N/A.

Variability in Global DNA Methylation Rate Across Tissues and Over Time in Sheep.

Recent studies highlighted the influence of epigenetic marks in the variability of many complex traits, both in plants and animals. These studied focused only on specific sites of the genome having differentially methylated profiles among individuals and/or tissues. In contrast, we recently used the methylation rate of the entire genome as a unique measure considered as a novel quantitative phenotype in sheep. This phenotype named global DNA methylation rate (GDMR), measured by luminometric assay, integrates the methylation level of each CpG dinucleotide within the 6 million of CCGG sites along the ovine genome. GDMR measured in blood previously showed moderate heritability of 0.20 and provided evidence for a genetic determinism. The main objective of the present study was to better characterize the GDMR phenotype in various tissues and investigate its variability in several breeds of sheep reared in the same environment. GDMR was measured on blood samples collected monthly from 59 growing male and female lambs (24 Romane, 23 Blackbelly and 12 Charollais), between birth and 4 months of age. Blood GDMR was on average around 80% and was influenced by the sampling date (p < 0.001), the breed (p = 0.002) and the sex (p = 0.002). In addition, GDMR was determined in 12 somatic (frontal lobe, pituitary gland, heart, lung, sub cutaneous and perirenal adipose tissue, skeletal muscle, liver, spleen, adrenal gland, medulla and cortical kidney) and 6 reproductive tissues (ovary, oviduct, uterus, testis, epididymis and seminal vesicle). GDMR was on average 70% in somatic tissues but marked variation was observed depending on the tissue. The GDMR measured in blood was higher than that measured in other somatic tissues, and is not a good proxy of less accessible tissues. Female reproductive tissues had a 10% higher GDMR than male reproductive tissues. We demonstrated a significant influence of the breed on blood GDMR, certainly reflecting the influence of different genetic backgrounds. The effect of the breed on GDMR may be related to their specific abilities to adapt to and live in different conditions.

Endothelial cell-derived fibroblast growth factor-18 regulates ovarian function in sheep.

Increasing the efficiency of farm animal reproduction is necessary to reduce the environmental impact of food production systems. One approach is to increase the number of healthy eggs (oocytes) produced per female for fertilization, thus it is important to understand factors that decrease oocyte health. One paracrine factor that decreases ovarian follicle growth is fibroblast growth factor 18 (FGF18) secreted by cells in the theca layer of the ovarian follicle, however the factors that regulate FGF18 secretion are unknown. In this study we hypothesized that FGF18 secretion is controled by intrafollicular factors and is linked to fertility, which we tested by using cell culture and sheep genetic models in vivo. Separation of theca cell populations revealed that FGF18 messenger RNA (mRNA) is located mainly in thecal endothelial rather than endocrine cells, and immunohistochemistry localized FGF18 protein to microvessels in the theca layer in situ. Culture of ovine theca-derived endothelial cells was used to demonstrate stimulation of FGF18 mRNA and protein abundance by bone morphogenetic protein 4 (BMP4), a growth factor derived from theca endocrine cells. Taking advantage of a sheep genetic model, we demonstrate reduced ovarian and peripheral FGF18 concentrations in the hyperprolific Booroola ewe harboring the FecBB mutation in BMPR1B. These data suggest a novel control of fertility by follicular endothelial cells, in which theca endocrine cells secrete BMP4 that stimulates the secretion of FGF18 from thecal endothelial cells, which in turn diffuses into the granulosa cell layer and promotes apoptosis.

Tissue Resources for the Functional Annotation of Animal Genomes.

In order to generate an atlas of the functional elements driving genome expression in domestic animals, the Functional Annotation of Animal Genome (FAANG) strategy was to sample many tissues from a few animals of different species, sexes, ages, and production stages. This article presents the collection of tissue samples for four species produced by two pilot projects, at INRAE (National Research Institute for Agriculture, Food and Environment) and the University of California, Davis. There were three mammals (cattle, goat, and pig) and one bird (chicken). It describes the metadata characterizing these reference sets (1) for animals with origin and selection history, physiological status, and environmental conditions; (2) for samples with collection site and tissue/cell processing; (3) for quality control; and (4) for storage and further distribution. Three sets are identified: set 1 comprises tissues for which collection can be standardized and for which representative aliquots can be easily distributed (liver, spleen, lung, heart, fat depot, skin, muscle, and peripheral blood mononuclear cells); set 2 comprises tissues requiring special protocols because of their cellular heterogeneity (brain, digestive tract, secretory organs, gonads and gametes, reproductive tract, immune tissues, cartilage); set 3 comprises specific cell preparations (immune cells, tracheal epithelial cells). Dedicated sampling protocols were established and uploaded in https://data.faang.org/protocol/samples. Specificities between mammals and chicken are described when relevant. A total of 73 different tissues or tissue sections were collected, and 21 are common to the four species. Having a common set of tissues will facilitate the transfer of knowledge within and between species and will contribute to decrease animal experimentation. Combining data on the same samples will facilitate data integration. Quality control was performed on some tissues with RNA extraction and RNA quality control. More than 5,000 samples have been stored with unique identifiers, and more than 4,000 were uploaded onto the Biosamples database, provided that standard ontologies were available to describe the sample. Many tissues have already been used to implement FAANG assays, with published results. All samples are available without restriction for further assays. The requesting procedure is described. Members of FAANG are encouraged to apply a range of molecular assays to characterize the functional status of collected samples and share their results, in line with the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles.

The ovarian follicle of ruminants: the path from conceptus to adult.

This review resulted from an international workshop and presents a consensus view of critical advances over the past decade in our understanding of follicle function in ruminants. The major concepts covered include: (1) the value of major genes; (2) the dynamics of fetal ovarian development and its sensitivity to nutritional and environmental influences; (3) the concept of an ovarian follicle reserve, aligned with the rise of anti-Müllerian hormone as a controller of ovarian processes; (4) renewed recognition of the diverse and important roles of theca cells; (5) the importance of follicular fluid as a microenvironment that determines oocyte quality; (6) the 'adipokinome' as a key concept linking metabolic inputs with follicle development; and (7) the contribution of follicle development to the success of conception. These concepts are important because, in sheep and cattle, ovulation rate is tightly regulated and, as the primary determinant of litter size, it is a major component of reproductive efficiency and therefore productivity. Nowadays, reproductive efficiency is also a target for improving the 'methane efficiency' of livestock enterprises, increasing the need to understand the processes of ovarian development and folliculogenesis, while avoiding detrimental trade-offs as greater performance is sought.

Identification of homozygous haplotypes carrying putative recessive lethal mutations that compromise fertility traits in French Lacaune dairy sheep.

BACKGROUND: Homozygous recessive deleterious mutations can cause embryo/fetal or neonatal lethality, or genetic defects that affect female fertility and animal welfare. In livestock populations under selection, the frequency of such lethal mutations may increase due to inbreeding, genetic drift, and/or the positive pleiotropic effects of heterozygous carriers on selected traits. RESULTS: By scanning the genome of 19,102 Lacaune sheep using 50 k single nucleotide polymorphism (SNP) phased genotypes and pedigree data, we identified 11 Lacaune deficient homozygous haplotypes (LDHH1 to LDHH11) showing a highly significant deficit of homozygous animals ranging from 79 to 100%. These haplotypes located on chromosomes 3, 4, 13, 17 and 18, spanned regions from 1.2 to 3.0 Mb long with a frequency of heterozygous carriers between 3.7 and 12.1%. When we compared at-risk matings (between carrier rams and daughters of carrier rams) and safe matings, seven of the 11 haplotypes were associated with a significant alteration of two fertility traits, a reduced success of artificial insemination (LDHH1, 2, 8 and 9), and/or an increased stillbirth rate (LDHH3, 6, 8, 9, and 10). The 11 haplotypes were also tested for a putative selective advantage of heterozygous carrier rams based on their daughter yield deviation for six dairy traits (milk, fat and protein yields, fat and protein contents and lactation somatic cell score). LDHH1, 3, 4, 5, 7, 9 and 11 were associated with positive effects on at least one selected dairy trait, in particular milk yield. For each haplotype, the most probable candidate genes were identified based on their roles in lethality of mouse knock-out models and in mammalian genetic disorders. CONCLUSIONS: Based on a reverse genetic strategy, we identified at least 11 haplotypes with homozygous deficiency segregating in French Lacaune dairy sheep. This strategy represents a first tool to limit at-risk matings in the Lacaune dairy selection scheme. We assume that most of the identified LDHH are in strong linkage disequilibrium with a recessive lethal mutation that affects embryonic or juvenile survival in sheep but is yet to be identified.

A Nonsense Variant in CCDC65 Gene Causes Respiratory Failure Associated with Increased Lamb Mortality in French Lacaune Dairy Sheep.

We recently demonstrated that the Lacaune deficient homozygous haplotype 6 (LDHH6) potentially hosts a recessive perinatal lethal mutation in Lacaune dairy sheep mapped on OAR3. In the present study, we have analyzed the whole-genome sequences of two Lacaune ram heterozygous carriers of LDHH6. After variant calling and filtering against the variants of 86 non-carrier rams, we have identified a single nucleotide variant (SNV) in the two LDHH6 carriers whose variant allele induced a premature stop codon (p.Glu111*) in the Coiled-Coil Domain Containing 65 (CCDC65) gene. CCDC65 is involved in the assembly of the nexin-dynein regulatory complex for the formation of microtubules in ciliated cells. In order to identify the phenotype in homozygous sheep, we generated at-risk matings (n = 17) between rams and ewes heterozygous for the candidate variant in CCDC65. A total of 16 lambs were born alive with five genotyped as homozygous carriers. The homozygous lambs suffered from respiratory problems, and four of them died within the first month of life. At necropsy, we observed a broad hepatization of lung lobes possibly induced by infectious pneumonia. The management of this lethal recessive allele (frequency of 0.06) through reasoned mating in the Lacaune sheep selection schemes could reduce lamb mortality by 2%.

New Anti-Müllerian Hormone Target Genes Involved in Granulosa Cell Survival in Women With Polycystic Ovary Syndrome.

PURPOSE: A protective effect of anti-Müllerian hormone (AMH) on follicle atresia was recently demonstrated using long-term treatments, but this effect has never been supported by mechanistic studies. This work aimed to gain an insight into the mechanism of action of AMH on follicle atresia and on how this could account for the increased follicle pool observed in women with polycystic ovary syndrome (PCOS). METHODS: In vivo and in vitro experiments were performed to study the effects of AMH on follicle atresia and on the proliferation and apoptosis of granulosa cells (GCs). RNA-sequencing was carried out to identify new AMH target genes in GCs. The expression of some of these genes in GCs from control and PCOS women was compared using microfluidic real time quantitative RT-PCR. RESULTS: A short-term AMH treatment prevented follicle atresia in prepubertal mice. Consistent with this result, AMH inhibited apoptosis and promoted proliferation of different models of GCs. Moreover, integrative biology analyses of 965 AMH target genes identified in 1 of these GC models, confirmed that AMH had initiated a gene expression program favoring cell survival and proliferation. Finally, on 43 genes selected among the most up- and down-regulated AMH targets, 8 were up-regulated in GCs isolated from PCOS women, of which 5 are involved in cell survival. MAIN CONCLUSIONS: Our results provide for the first time cellular and molecular evidence that AMH protects follicles from atresia by controlling GC survival and suggest that AMH could participate in the increased follicle pool of PCOS patients.

Genome-Wide Identification of a Regulatory Mutation in BMP15 Controlling Prolificacy in Sheep.

The search for the genetic determinism of prolificacy variability in sheep has evidenced several major mutations in genes playing a crucial role in the control of ovulation rate. In the Noire du Velay (NV) sheep population, a recent genetic study has evidenced the segregation of such a mutation named FecL L . However, based on litter size (LS) records of FecL L non-carrier ewes, the segregation of a second prolificacy major mutation was suspected in this population. In order to identify this mutation, we have combined a case/control genome-wide association study with ovine 50k SNP chip genotyping, whole genome sequencing, and functional analyses. A new single nucleotide polymorphism (OARX:50977717T > A, NC_019484) located on the X chromosome upstream of the BMP15 gene was evidenced to be highly associated with the prolificacy variability (P = 1.93E-11). The variant allele was called FecX N and shown to segregate also in the Blanche du Massif Central (BMC) sheep population. In both NV and BMC, the FecX N allele frequency was estimated close to 0.10, and its effect on LS was estimated at +0.20 lamb per lambing at the heterozygous state. Homozygous FecX N carrier ewes were fertile with increased prolificacy in contrast to numerous mutations affecting BMP15. At the molecular level, FecX N was shown to decrease BMP15 promoter activity and supposed to impact BMP15 expression in the oocyte. This regulatory action was proposed as the causal mechanism for the FecX N mutation to control ovulation rate and prolificacy in sheep.

A new allele in the BMP15 gene (FecXRA) that affects prolificacy co-segregates with FecXR and FecXGR in Rasa aragonesa sheep.

A FecX-mutated allele called FecXR in the BMP15 gene has been described in Rasa aragonesa sheep. FecXR causes increased prolificacy when heterozygous and sterility when homozygous in ewes. However, highly prolific ewes without the FecXR allele have been found in this breed. Therefore, a genome-wide association study (GWAS) was performed in 158 ewes (tail H: N = 73, mean prolificacy ± standard deviation = 2.14 ± 0.26; tail L: N = 85, mean prolificacy = 1.06 ± 0.08) with the Ovine HD SNP BeadChip. In this analysis, the FecXGR allele was found to have genome-wide significance associated with prolificacy, first described in the Grivette sheep breed. We also identified a novel polymorphism in exon 2 of BMP15 in 9 high prolific ewes by Sanger sequencing. This new mutation, called FecXRA, is a SNP (Oar3.1_X: g. 50970948C > T; NM_001114767.1: c.1172C > T) that produces an amino acid substitution (ENSOART00000010201: p.T400I) that is predicted to be deleterious and to alter the predicted secondary structure of the mature protein. To confirm if this SNP had any the effect on prolificacy, we genotyped sires with known EBVs (Estimated Breeding Values), finding one hemizygous sire for the FecXRA allele with the highest EBV in the breeding program (effect on litter size at + 0.39 lamb per lambing). A very low frequency, ranging from 0.13 to 2%, was found for the FecXGR and FecXRA alleles in 3428 animals belonging to four different flocks. Finally, an association study was performed to validate and quantify the effects of the FecXGR and FecXRA alleles. Significant increased prolificacy of 0.52 ± 0.05, 0.42 ± 0.05 and 0.32 ± 0.01 were found when comparing FecXGR, FecXRA and FecXR heterozygous ewes to wild type homozygous ones. These effects are of the same order of magnitude as the effect of most of other known major genes for prolificacy. Only significant differences between FecXGR and FecXR were found among the three alleles associated with increased prolificacy. However, we cannot confirm the effect of the FecXRA allele at homozygous state because we did not find any homozygous ewes. These results confirm that these three alleles in the BMP15 gene that affect prolificacy co-segregate in Rasa aragonesa sheep.

Prenatal programming by testosterone of follicular theca cell functions in ovary.

In mammalian ovaries, the theca layers of growing follicles are critical for maintaining their structural integrity and supporting androgen synthesis. Through combining the postnatal monitoring of ovaries by abdominal magnetic resonance imaging, endocrine profiling, hormonal analysis of the follicular fluid of growing follicles, and transcriptomic analysis of follicular theca cells, we provide evidence that the exposure of ovine fetuses to testosterone excess activates postnatal follicular growth and strongly affects the functions of follicular theca in adulthood. Prenatal exposure to testosterone impaired androgen synthesis in the small antral follicles of adults and affected the expression in their theca cells of a wide array of genes encoding extracellular matrix components, their membrane receptors, and signaling pathways. Most expression changes were uncorrelated with the concentrations of gonadotropins, steroids, and anti-Müllerian hormone in the recent hormonal environment of theca cells, suggesting that these changes rather result from the long-term developmental effects of testosterone on theca cell precursors in fetal ovaries. Disruptions of the extracellular matrix structure and signaling in the follicular theca and ovarian cortex can explain the acceleration of follicle growth through altering the stiffness of ovarian tissue. We propose that these mechanisms participate in the etiology of the polycystic ovarian syndrome, a major reproductive pathology in woman.

Genetic Determinism Exists for the Global DNA Methylation Rate in Sheep.

Recent studies showed that epigenetic marks, including DNA methylation, influence production and adaptive traits in plants and animals. So far, most studies dealing with genetics and epigenetics considered DNA methylation sites independently. However, the genetic basis of the global DNA methylation rate (GDMR) remains unknown. The main objective of the present study was to investigate genetic determinism of GDMR in sheep. The experiment was conducted on 1,047 Romane sheep allocated into 10 half-sib families. After weaning, all the lambs were phenotyped for global GDMR in blood as well as for production and adaptive traits. GDMR was measured by LUminometric Methylation Analysis (LUMA) using a pyrosequencing approach. Association analyses were conducted on some of the lambs (n = 775) genotyped by using the Illumina OvineSNP50 BeadChip. Blood GDMR varied among the animals (average 70.7 ± 6.0%). Female lambs had significantly higher GDMR than male lambs. Inter-individual variability of blood GDMR had an additive genetic component and heritability was moderate (h 2 = 0.20 ± 0.05). No significant genetic correlation was found between GDMR and growth or carcass traits, birthcoat, or social behaviors. Association analyses revealed 28 QTLs associated with blood GDMR. Seven genomic regions on chromosomes 1, 5, 11, 17, 24, and 26 were of most interest due to either high significant associations with GDMR or to the relevance of genes located close to the QTLs. QTL effects were moderate. Genomic regions associated with GDMR harbored several genes not yet described as being involved in DNA methylation, but some are already known to play an active role in gene expression. In addition, some candidate genes, CHD1, NCO3A, KDM8, KAT7, and KAT6A have previously been described to be involved in epigenetic modifications. In conclusion, the results of the present study indicate that blood GDMR in domestic sheep is under polygenic influence and provide new insights into DNA methylation genetic determinism.

Multi-species annotation of transcriptome and chromatin structure in domesticated animals.

BACKGROUND: Comparative genomics studies are central in identifying the coding and non-coding elements associated with complex traits, and the functional annotation of genomes is a critical step to decipher the genotype-to-phenotype relationships in livestock animals. As part of the Functional Annotation of Animal Genomes (FAANG) action, the FR-AgENCODE project aimed to create reference functional maps of domesticated animals by profiling the landscape of transcription (RNA-seq), chromatin accessibility (ATAC-seq) and conformation (Hi-C) in species representing ruminants (cattle, goat), monogastrics (pig) and birds (chicken), using three target samples related to metabolism (liver) and immunity (CD4+ and CD8+ T cells). RESULTS: RNA-seq assays considerably extended the available catalog of annotated transcripts and identified differentially expressed genes with unknown function, including new syntenic lncRNAs. ATAC-seq highlighted an enrichment for transcription factor binding sites in differentially accessible regions of the chromatin. Comparative analyses revealed a core set of conserved regulatory regions across species. Topologically associating domains (TADs) and epigenetic A/B compartments annotated from Hi-C data were consistent with RNA-seq and ATAC-seq data. Multi-species comparisons showed that conserved TAD boundaries had stronger insulation properties than species-specific ones and that the genomic distribution of orthologous genes in A/B compartments was significantly conserved across species. CONCLUSIONS: We report the first multi-species and multi-assay genome annotation results obtained by a FAANG project. Beyond the generation of reference annotations and the confirmation of previous findings on model animals, the integrative analysis of data from multiple assays and species sheds a new light on the multi-scale selective pressure shaping genome organization from birds to mammals. Overall, these results emphasize the value of FAANG for research on domesticated animals and reinforces the importance of future meta-analyses of the reference datasets being generated by this community on different species.

The high prolificacy of D'man sheep is associated with the segregation of the FecLL mutation in the B4GALNT2 gene.

Mutations in the FecL locus are associated with large variation in ovulation rate and litter size in the French Lacaune sheep breed. It has been shown that the B4GALNT2 gene within the FecL locus is most likely responsible for the high fecundity in the French breed. In this study, we have highlighted the segregation of the FecLL mutation within the B4GALNT2 gene in North African sheep breeds and notably in the highly prolific D'man breed. Genotyping of a sample of 183 Tunisian D'man individuals revealed a high frequency (0.65) of the prolific allele FecLL which was attributed to the adoption of a decades-old breeding strategy based on the selection of ewe lambs born from large litter size. Homozygous LL ewes showed a significantly increased litter size compared to heterozygous and non-carrier ewes (FecLL /FecLL  = 2.47 ± 0.09 vs. FecLL /FecL+  = 2.23 ± 0.09, p < 0.05 and FecL+ /FecL+  = 1.93 ± 0.18, p < 0.01). The presence of the FecLL polymorphism in both D'man and Lacaune breeds argues for an ancestral origin of this mutation and brings an answer to the old question of the genetic determinism of the extreme prolificacy of the D'man ewes. The results of this study can help to establish planned genotype-based mating allowing both higher profit for the breeders and an optimal management of the FecLL mutation in D'man sheep populations.

Transcriptome analysis of ovine granulosa cells reveals differences between small antral follicles collected during the follicular and luteal phases.

Ovarian follicular growth occurs in both the follicular and luteal phases of the estrous cycle but in very different endocrine contexts. In both phases, many small antral follicles with similar morphologic and histologic characteristics are present within the ovaries as a reserve for the terminal folliculogenesis. However, there are several gaps in our molecular knowledge of the gene expression profiles of small antral follicles in the follicular and luteal phases. The aim of the present study was to use RNA sequencing to compare and analyze the global transcriptional profile of ovine granulosa cells collected from small antral follicles (1-3 mm) either during the follicular or the luteal phase of the estrous cycle, with the hypothesis that they should be differential. We identified 663 genes whose mRNA was differentially expressed or accumulated in the granulosa cell layer of small antral follicles in the two phases. A comprehensive interpretation of these data was performed through integrative analyses (Gene Ontology, Ingenuity Pathway Analysis) and the exploitation of already available transcriptomic data on follicular growth and atresia. In particular, we observed that the contrasted endocrine context between follicular and luteal phases may have an impact on estradiol, follicle-stimulating hormone (FSH), and on the activin/inhibin signaling pathways. Furthermore, we reveal the possible initiation of early follicular atresia in small antral follicles during the follicular phase in interaction with the presence of immune cells. This study provides new insights into the gene expression profile in ovine granulosa cells, and we suggest that these molecular changes may have an implication at the time of follicle selection.

Variation in Recombination Rate and Its Genetic Determinism in Sheep Populations.

Recombination is a complex biological process that results from a cascade of multiple events during meiosis. Understanding the genetic determinism of recombination can help to understand if and how these events are interacting. To tackle this question, we studied the patterns of recombination in sheep, using multiple approaches and data sets. We constructed male recombination maps in a dairy breed from the south of France (the Lacaune breed) at a fine scale by combining meiotic recombination rates from a large pedigree genotyped with a 50K SNP array and historical recombination rates from a sample of unrelated individuals genotyped with a 600K SNP array. This analysis revealed recombination patterns in sheep similar to other mammals but also genome regions that have likely been affected by directional and diversifying selection. We estimated the average recombination rate of Lacaune sheep at 1.5 cM/Mb, identified ∼50,000 crossover hotspots on the genome, and found a high correlation between historical and meiotic recombination rate estimates. A genome-wide association study revealed two major loci affecting interindividual variation in recombination rate in Lacaune, including the RNF212 and HEI10 genes and possibly two other loci of smaller effects including the KCNJ15 and FSHR genes. The comparison of these new results to those obtained previously in a distantly related population of domestic sheep (the Soay) revealed that Soay and Lacaune males have a very similar distribution of recombination along the genome. The two data sets were thus combined to create more precise male meiotic recombination maps in Sheep. However, despite their similar recombination maps, Soay and Lacaune males were found to exhibit different heritabilities and QTL effects for interindividual variation in genome-wide recombination rates. This highlights the robustness of recombination patterns to underlying variation in their genetic determinism.